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Our low-temperature magneto-Raman scattering measurements performed on graphene-like loca¬ 
tions on the surface of bulk graphite reveal a new series of magneto-phonon resonances involving 
both K-point and F point phonons. In particular, we observe for the first time the resonant splitting 
of three crossing excitation branches. We give a detailed theoretical analysis of these new reso¬ 
nances. Our results highlight the role of combined excitations and the importance of multi-phonon 
processes (from both K and F points) for the relaxation of hot carriers in graphene. 

PACS numbers: 73.22.Lp, 63.20.Kd, 78.30.Na, 78.67.-n 


I. INTRODUCTION 

Quantization of electronic motion by a strong magnetic 
field leads to many qualitatively new effects. They are 
most pronounced in two-dimensional systems, where the 
motion along the field is frozen by quantum confinement, 
so the continuous electronic spectrum is transformed into 
a series of discrete Landau levels (LLs). Probably the 
most famous examples of such new effects are the inte¬ 
ger and the fractional quantum Hall effects [T], which 
arise because the degenerate electronic states on a LL 
are profoundly modified by disorder or lateral confine¬ 
ment potential, and by electron-electron interaction, re¬ 
spectively. Electron-phonon coupling can also produce 
new effects, especially when the energy of an inter-LL 
excitation matches that of an optical phonon. Such res¬ 
onances may lead to strong-coupling effects even if the 
electron-phonon coupling in the absence of the magnetic 
field was weak. 

Magneto-phonon resonances were first discussed theo¬ 
retically in relation to the dc electron transport in doped 
bulk semiconductors subject to a strong magnetic field 
where the resonance leads to an enhancement of the 
electron-phonon scattering rate [5]. They were later ob¬ 
served experimentally in InSb Em]. A doublet struc¬ 
ture was observed in the infrared absorption spectrum of 
bulk InSb in a magnetic field, which was attributed to 
the resonant coupling of an electronic excitation, the cy¬ 
clotron resonance mode, to optical phonons, when tuned 
in resonance [5]. This is a strong-coupling effect, whose 
detailed theory was presented in Ref. [5]. Later, it was 
suggested that such magneto-phonon resonance need not 
be restricted to doublets: three branches of excitations 
involving zero, one, and two phonons may also cross at 
the same value of the magnetic field, resulting in a triplet 
structure of the absorption peak [7]. In fact, equidistancy 
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of LLs for electrons with parabolic bands leads to a rich 
variety of possible resonant combinations [8]. Still, we 
are not aware of any observation of a resonant splitting of 
more than two excitation branches in conventional semi¬ 
conductors. 

Graphene, with its conical electron dispersion and ex¬ 
ceptional crystal quality, opens new possibilities for ob¬ 
servation of fundamental physical effects. In a magnetic 
field B, perpendicular to the crystal plane, the LL ener¬ 
gies are E±n = ±V^{^v/Ib), where n = 0,1,... is the 
Landau level index, v « 10®m/s is the Dirac velocity, and 
Ib = \/fi/(eB) is the magnetic length. Due to the y/n de¬ 
pendence, the spacing between LLs with small n is larger 
than for electrons with parabolic spectrum, which en¬ 
abled the observation of the quantum Hall effect at room 
temperature [9], unprecedented in conventional semicon¬ 
ductor structures. Magnetophonon resonances have been 
observed in Raman scattering on optical phonons mm, 
where the Raman G peak acquired a doublet structure. 
Here, we report an observation of magnetophonon reso¬ 
nances where the triple peak structure is clearly seen. To 
the best of our knowledge, it is the first observation of 
such resonances in solid-state physics. 

Despite the generality of the observed magnetophonon 
resonance effects, graphene introduces some qualitatively 
new aspects. First, electronic bands in graphene are con¬ 
ical rather than parabolic, so the LLs are not equally 
spaced in energy. This leads to a classification of possi¬ 
ble resonant combinations which is quite different from 
that in conventional semiconductors. Such classification 
for graphene is developed below. Second, while in an 
undoped semiconductor all LLs in the conduction band 
are empty, and those in the valence band are completely 
filled, in graphene there is always a partially filled LL (un¬ 
less the magnetic field is tuned to some special values); if 
this partially filled level is involved in the resonance, the 
physics of the electron-phonon coupling is changed sig¬ 
nificantly. Namely, it turns out that the multiple-peak 
structure cannot be viewed as a resonant splitting of a 
few discrete levels; one has to face a true many-body 
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problem, so the peak splitting is necessarily accompa¬ 
nied by a broadening of a similar magnitude, and compli¬ 
cated spectral shapes may be produced. This physics was 
overlooked in the previous theoretical studies of magneto¬ 
phonon resonances in graphene [IMl]- 

Our experiments are performed on extremely pure 
flakes of graphene which can be found on the surface of 
bulk graphite [nunHH!. To probe the magnetophonon 
resonances, we use Raman scattering, which is a powerful 
and popular tool for characterization of graphene PO] , 
Inter-LL electronic excitations can be observed in the 
Raman spectrum of graphene [HI El], and the strongest 
features are due to transitions —n —> n [22j . In addition, 
transitions —n —)■ n±l are also observed [15] even though 
theory predicts them to be weaker |22j . In this work, we 
trace the evolution of inter-LL excitations in graphene 
on graphite in magnetic fields up to R = 30 T. Such 
fields are needed to tune the inter-LL excitation energies 
across the energies of optical phonons and their combi¬ 
nations. We observe different series of magnetophonon 
resonances, involving both phonons from the vicinity of 
the K and T points of the first Brillouin zone. 

The paper is organized as follows. In Sec. [IT] we dis¬ 
cuss qualitatively different types of magnetophonon reso¬ 
nances which can occur in graphene, and summarize our 
main results. In Sec. |III| we give a detailed presentation 
of our experimental results. In Sec. |IV| the theory is de¬ 
veloped for the various cases discussed in Sec. |ll] 


II. QUALITATIVE DISCUSSION: 

CLASSIFICATION OF MAGNETO-PHONON 
RESONANCES AND SUMMARY OF MAIN 
RESULTS 

In the following, we use a shorthand notation L-n,m 
for an excitation of an electron from a filled level —n 
to an empty level m (assuming that the partially filled 
level is n = 0). To gain a qualitative insight into various 
resonant processes, it is instructive to analyze the weak 
coupling limit first. Namely, consider an initial excitation 
which was created in the sample in the course of Raman 
scattering, and analyze different channels of decay for 
this excitation (see Fig.[^. 

One case is when the initial excitation is a phonon; 
for example, the doubly degenerate E 2 g phonon with the 
wave vector at the T point of the first Brillouin zone, giv¬ 
ing rise to the G peak in the Raman spectrum. It may 
decay into an electronic inter-LL excitation, as schemat¬ 
ically shown in Fig. [^a). The selection rule for the ma¬ 
trix element is n —> n ± 1, so the phonon decay rate is 
strongly enhanced when the phonon is resonant with one 
of the corresponding transitions, that is, the phonon en¬ 
ergy hwT = En+i - E^n = En - for soiue n. 

This case corresponds to the magneto-phonon resonance 
in the Raman scattering on phonons, which was observed 
experimentally mm, following the theoretical predic¬ 
tion in Refs. mm- 
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Figure 1. (color online) Different types of magneto-phonon 
resonance in excitation decay, (a) An optical phonon decay¬ 
ing into an electronic excitation, (b) An electronic excitation 
decaying into an optical phonon and another electronic exci¬ 
tation. The horizontal lines represent the Landau levels, the 
black circles represent their filling. The filled red circle repre¬ 
sents the excited electron, the empty red circle represents the 
hole, and the wavy line represents the phonon. 
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In the present work, we focus on a different case, that of 
the initial excitation being an inter-LL excitation which 
decays into a phonon and another inter-LL excitation. 
Again, the corresponding decay rate is enhanced at res¬ 
onance, but the resonance condition is different from the 
previous case. For example, if the initial excitation is 
L-n,n-, it can decay if the electron on the level n or the 
hole on the level —n emits a phonon and moves to an¬ 
other level n' or —n', respectively (n' < n). The de¬ 
cay is resonant if the phonon energy Huj = En — En'. 
It is essentially the same process that is responsible for 
the magneto-phonon resonance in electronic transport in 
semiconductors US], as well as in magneto-optical ab¬ 
sorption mi. 

The role of momentum conservation is quite different in 
the two processes. As we focus here on optical excitation, 
the momentum of the incident and scattered photons can 
be neglected, so in both cases, the total momentum of the 
initial excitation is zero. If the initial excitation is a single 
phonon, it can be only the E 2 g optical phonon from the 
F point, and the total momentum of the inter-LL excita¬ 
tion into which it decays should also be zero, which leads 
to the selection rule n —)■ nil, and thus restricts this case 
to the magneto-phonon resonance studied in Refs. nDi- 
HSj. For a multi-phonon initial excitation, one can also 
think about the resonant enhancement of the phonon de¬ 
cay rate; however, the width of the multi-phonon Raman 
peaks is dominated by other factors than the phonon de¬ 
cay [23U25] , so the modification of latter by the magnetic 
field is not observed in the Raman spectra. 

For an initial electronic inter-LL excitation of zero total 
momentum, the momentum conservation allows the de- 
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cay into a phonon with an arbitrary momentum and an¬ 
other inter-LL excitation with the opposite momentum. 
Still, the matrix element is strongly suppressed unless 
the phonon momentum is within ^1/1^ from the F or 
K, K' points of the first Brillouin zone. Thus, one has to 
consider two kinds of resonances: (i) those involving the 
E 2 g optical phonons around the F point and intravalley 
electronic inter-LL excitations, and (ii) those involving 
intervalley inter-LL excitations and phonons around the 
K, K' points. Among the latter phonons, those most 
strongly coupled to the electrons are the Aig transverse 
optical phonons. The corresponding phonon en ergi es are 
hjjJT — 196 meV and hujK = 160 meV (see Sec. III). An¬ 
other important consequence of finite phonon momenta 
1 /Ib) involved in the process, is the absence of any 
selection rule on the LL indices. In the above example of 
the decay of the L_n,n excitation into a phonon and the 
L-n,n', any n' is allowed (energy conservation requires 


n' < n). 

The above discussion of decay of the initial excitation 
is valid in the limit when the electron-phonon coupling is 
weak compared to the typical broadening of the Landau 
levels or to the dispersion of optical phonons or inter- 
LL excitations on the momentum scale 1/Ib- In the 
opposite limiting case, the problem corresponds to the 
coupling between discrete energy levels, when the states 
become strongly mixed. The energies of these hybrid 
excitations (magnetopolarons) exhibit an anticrossing as 
a function of the magnetic field as the latter is swept 
through the resonance region. 

For an initial electronic inter-LL excitation, different 
types of resonances are possible. A resonance occurs 
when the phonon energy hui matches the energy differ¬ 
ence between some empty Landau levels: 


^n' — 5 


( 1 ) 


for some n > n' ^ 0 (we assume to be at zero tem¬ 
perature and zero doping, so n = 0 is the only partially 
filled level). Due to electron-hole symmetry, = —En, 
Eq. 0 automatically yields a similar condition for two 
filled levels, E_n' — E_n = huj. Also, because of the -y/n 
dependence of En, the LLs are not equidistant, but still, 
some additional degeneracies may arise. For example, if 
E 4 — El = tvjj, this automatically yields Ei — Eq = hui, 
while ii Eq — E 2 = fvjj, there are no additional degenera¬ 
cies. 

All resonances of the kind 0 can be classified into 
three groups, depending on whether the ratio 


^ _ En _ _ 

Huj y/n - 


( 2 ) 


is integer, rational, or irrational (see Fig. 0. Let us see 
what happens in each of these cases, restricting our at¬ 
tention to the initial excitations of the type which 

dominate the electronic Raman scattering in the mag¬ 
netic field [22], and to L_(^n+i),n (as well as its electron- 
hole-symmetric counterpart which are also ob¬ 

served in the Raman spectra, as discussed in Sec. Ini 
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Figure 2. (color online) Examples of different conse¬ 
quences of a resonance of the type of Eq. 0. (i) Irrational 
^ — y/n/(y/n — 2) for the resonance En = E 4 -\- hui. No 
further resonant steps are possible, (ii) Rational (^ = 5/2 for 
the resonance E 25 = Eq + hui, which automatically yields the 
next resonant step, E 25 = Ei -\- 2huj. The n = 0 level is not 
involved, (iii) Integer C = 2 for the resonance Eig = E 4 -f hui, 
which yields the next resonant step, Eig = Eq + 2 huj, in¬ 
volving the partially filled zero level. The horizontal lines 
represent the Landau levels, the black circles represent their 
filling. The filled red circles represent the excited electrons, 
the vertical arrows show their transitions which are resonant 
with a phonon. 


(i) If C is irrational, the resonances do not proliferate. 
Then the simplest case is that of the initial 
excitation which is resonant with L_(^n+i),n' plus one 
phonon. In this case, there are no hole resonances; in¬ 
deed, E-n" — E_(^n+i) = Hui aX the same time would 
imply y/rT+A — y/n = y/rE — yjn' for some n! < n" < n, 
which is impossible due to the concavity of the square 
root. The symmetric excitation L_n,n is then resonant 
with L-n,n' (or L-n',n) plus One phonon, but also with 
L-n’,n' plus two phonons, i. e. the electron and the hole 
resonances necessarily occur simultaneously. As the ex- 
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citations with zero, one, and two phonons disperse differ¬ 
ently with the magnetic field, the symmetric excitation 
gives a triple anticrossing of the levels as a function of the 
magnetic field. Such triple anticrossings were discussed 
long ago in the context of optical magnetoabsorption in 
bulk semiconductors [7]. 

(ii) Let C be rational but not an integer, and let m be 

the integer part of C- In this case, condition 0 yields 
not just one, but m resonances. Indeed, for any j ^ m 
the energy En — jhui matches the energy En^ of the level 
rij = [j\/W — {j — The last level in this se¬ 

quence has the energy En^ = (C ~ mjSw, which does 
not allow for emission of any more phonons. However, 
if we construct n^-i-i = [{m + l)^/n' — then it 

turns out that + y/rim+i = y/n — or, equiva¬ 
lently, huj = En^ — E-nm+i ■ This means that the phonon 
itself is resonantly hybridized with the electronic excita¬ 
tion from the filled level —rim+i to the empty level rim 
(as the phonon involved has a finite wave vector, there 
are no selection rules on nm,rim+i)- Thus, the asym¬ 
metric excitation L_(n±i),n is resonant with 

plus j phonons for any j = 1,..., m, so we have an an¬ 
ticrossing of TO -|- 1 energy levels with different disper¬ 
sions (the hybridization of the phonon with 
or L-n^,n^^i does not add anything, as these electronic 
excitations disperse with the magnetic field in the same 
way as En — En'). The symmetric excitation L-n,n can 
resonantly emit phonons both on the electron and hole 
sides, up to 2 to phonons in total, corresponding to the 
anticrossing of 2m + 1 levels. It is important that for 
non-integer ^ none of the resonances involve the partially 
filled n = 0 Landau level. 

(iii) If C = m, an integer, then the last term of the 
sequence {nj} from the previous case is rim = 0. The 
resonance counting is analogous to case (ii). Namely, we 
have TO -|- 1 anticrossings for the asymmetric initial ex¬ 
citation L_(„±i),„ and 2 m + 1 anticrossings for 

so that the latter can thus be entirely transformed into 
2m phonons by a sequence of resonant transitions. How¬ 
ever, involvement of the partially filled n = 0 level brings 
about a qualitatively new feature: the possibility of cre¬ 
ation of an arbitrary number of zero-energy electronic ex¬ 
citations on the partially filled level. As will be discussed 
in Sec. HYd] this feature makes a consistent theoretical 
treatment of case (iii) quite problematic. 

Our experimental results, presented in Sec. [ml include 
a clear observation of several resonances of type (iii) with 
double and, for the first time, triple avoided crossings. 
Near each avoided crossing, we describe our data phe¬ 
nomenologically by an effective two- or three-level model 
with some coupling between the resonant excitations. We 
also observe some signatures of resonances of type (i), 
which are, however, close to the limits of our experimen¬ 
tal resolution, and they do not allow for a quantitative 
anlysis. Resonances of type (ii) are quite difficult to ob¬ 
serve, as the smallest Landau level indices n, n' which 
produce non-integer rational C are n = 9, n' = 1, leading 
to C = 3/2. The corresponding magnetic field is too low. 


and the Landau levels involved in the electronic excita¬ 
tion are too high, to be resolved in our experiment. 

The theory for all three cases is presented in Sec. EYi 
Its main task is to justify (or disprove) the phenomeno¬ 
logical description of each avoided crossing by an effective 
few-level model. Clearly, such description may be correct 
only if one can neglect the energy dispersion of the elec¬ 
tronic excitations and of the phonons. The former can be 
neglected if the Coulomb interaction is sufficiently weak 
(e.g., due to screening by the conducting graphite sub¬ 
strate). As for the latter, one should recall that the most 
important contribution to the dispersion of phonons near 
the r or AT points comes from their coupling to the Dirac 
electrons which produces the Kohn anomaly [13 HZ]- In 
a strong magnetic field the spectrum of Dirac electrons 
is quantized, so the main effect of the electron-phonon 
coupling is to induce magneto-phonon resonances. The 
residual part of the phonon dispersion which is due to 
interaction with electrons in the a bands is quite weak 
and can indeed be neglected, as the relevant scale of the 
phonon wave vectors is quite small, q ^ I/Ib {^b is the 
magnetic length). 

Still, it turns out, that even if one starts from dis¬ 
persionless excitations, energy dispersion can be eventu¬ 
ally generated due to resonant coupling. This effect is 
especially dramatic for case (iii), where plenty of zero- 
energy excitations are available on the partially filled 
n = 0 Landau level. As a result, the splitting between the 
energies of the coupled excitations is accompanied by a 
broadening, whose magnitude is of the same order as the 
splitting. The broadened peaks have complicated spec¬ 
tral shapes, not described by any simple functional form. 


such as Lorentzian, Gaussian, or any other (see Fig. 11 


Sec. IVD for an example). Thus, the phenomenological 
description of our experimental data by effective two- or 
three-level models is very approximate. 

The splitting of the Lo,i electronic excitation, observed 
when it is in resonance with the T point phonon near 
H = 28 T [ini HI] , can be ascribed to either of the two 
types of resonances mentioned in the beginning of this 
section. Indeed, on the one hand, it can be viewed as the 
decay of the initial phonon into an electronic excitation. 
On the other hand, it can also be viewed as the decay of 
the initial Lq.i electronic excitation into a phonon and a 
zero-energy Lq.o excitation, falling into case (iii) above. 
Due to such “double-faced” nature of this particular res¬ 
onance, the split peaks are strongly broadened. This 
intrinsic broadening is missing in the simple description 
of Refs. [T2l[l3]. For the resonance between the AT point 
phonon and the Lq i excitation, coupling to an exter¬ 
nal continuum of states was introduced in Ref. [55| to 
describe the strong broadening of the peaks; here we see 
that the broadening, in fact, originates from the electron- 
phonon coupling itself, and there is no need to invoke any 
extrinsic effects. 


The broadening effects are less dramatic in the cases 
(i) and (ii). As discussed above, case (ii) necessarily 
involves Landau levels with quite large n ^ 9. This 
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strongly suppresses the broadening effect: the main peaks 
are narrow and they contain the most of the spectral 
weight (Sec. IVC I. In case (i), discussed in Sec. IVB for 
asymmetric transitions Ij_(n±i),n no broadening is gen¬ 
erated at all, so their resonance with i_(n±i),n' and a 
phonon can be faithfully described by an effective two- 
level model. No broadening arises also for symmetric 
transitions L_n,n resonating with and the K point 

phonon, so the effective three-level description is justi¬ 
fied. However, if F point phonons are involved in the res¬ 
onance, broadening is generated due to a peculiar phonon 
exchange effect: the phonons emitted by i-n,o and Lo,n 
excitations which constitute the initial F_„,n excitation, 
can be reabsorbed in the opposite order; this is impossi¬ 
ble for the K point phonons as the two emitted phonons 
must necessarily belong to the opposite valleys. 


III. EXPERIMENT 


The low-temperature magneto-Raman scattering re¬ 
sponse, at the pm scale, has been measured using a home¬ 
made miniaturized optical probe based on optical fibers 
(a 5 pm core mono-mode fiber for the excitation and a 
50 pm core for the collection), lenses and optical band 
pass filters to clean the laser and to reject the laser line. 
The excitation laser at A = 514.5 nm is focused on the 
sample with a high numerical aperture aspherical lens 
down to a spot of the diameter ^ 1 pm. The sample 
is placed on piezo stages to allow for the spatial map¬ 
ping of the Raman response of the sample. This optical 
set-up is inserted in a closed jacket, with some helium 
exchange gas, which is then immersed in liquid helium 
and placed in the center of resistive solenoid producing 
magnetic fields up to 30 T. The unpolarized Raman scat¬ 
tering response of our samples has been measured in the 
quasi-backscattering geometry with the magnetic field B 
applied perpendicularly to the graphene crystal plane. 
This experimental set-up has been described in more de¬ 
tails in Refs. [laiin]. To find a graphene flake on the 
surface of bulk graphite, we use the methodology based 
on the spatial mapping of the Raman scattering response 
of the surface of bulk graphite with an applied magnetic 
field described in Ref. [HO] . 

The typical Raman spectra, from which the i? = 0 re¬ 
sponse has been subtracted, are shown in Fig. [^a). The 
two gray bars mask the two strong phonon features, the 
G band around 1582 cm“^ and the 2D band observed at 
2750 cm“^ for our 514.5 nm excitation. Typically, the ob¬ 
served line widths for electronic features are ~ 30 cm“^. 
Due to the pronounced difference in the scattered in¬ 
tensity of phonon excitations and of the different types 
of electronic excitations, we present our main results in 
Fig. §b) in the form of a false color map of the scat¬ 
tered intensity as a function of the magnetic field which 
has been differentiated with respect to the magnetic field, 
with a step of 100 mT. The 2D band feature is weakly 
affected by the magnetic field, but because of a Faraday 
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Figure 3. a) Characteristic Raman scattering spectra at 
selected values of the magnetic field, from which the zero- 
field response has been subtracted. The vertical gray bars 
mask the residual, after normalization, phonon contributions, 
due to the Faraday effect in the experimental set-up. b) False 
color map of the Raman scattering response of the graphene¬ 
like location, differentiated with respect to B, as a function of 
the magnetic field (with the step in B of 100 mT). Main L-n,m 
electronic excitations are indicated. The horizontal solid line 
represents Twjk, the dashed line 2hLOK, the dashed-dotted line 
is at the energy position of the 2D band feature while the 
dotted line represents 2hoJr- The three white (black) boxes 
identify spectral and magnetic field regions around B — 10,15 
and 28T (B = 6, 9 and 19T) where the electron-F point (K 
point) phonon resonance affects the spectra. 


effect affecting the intensity of the excitation laser when 
sweeping the magnetic field, it cannot be cleanly removed 
by normalizing spectra acquired at a finite field by the 
B = 0 response. Thus, differentiation as a function of B 
provides the cleanest color map. 

The rich series of avoided crossings around 1582 cm“^ 
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magnetic field 

resonant levels 

transitions affected 

28.5 T 

El = ficur 

±0.1, ±-1,1, ±-1,2 

18.7 T 

E\ = huJK 

±0,1, ±-1,1, ±-1,2 

14.2 T 

E2 — hLOY 

±-1,2, ±-2,2, ±-2,3 

9.5 T 

3 

II 

±-2,3, ±-3,3, ±-3,4 

9.4 T 

±2 = hijJK 

±-1,2, ±-2,2, ±-2,3 

6.2 T 

±3 = hLOK 

±-2,3, ±-3,3, ±-3,4 


Table I. The values of the magnetic fields, evaluated theoreti¬ 
cally, and the corresponding resonances which can be observed 


in Fig.^b). 


is due to the well-studied magneto-phonon resonance in 
the Raman scattering on the F-point optical phonon |10F 
m- In the following, instead of tracing the evolution of 
the phonon feature with B, we study the evolution of 
the inter-LL electronic excitations. This approach is to¬ 
day only possible for graphene on graphite which shows 
an exceptionally high quality [TB], surpassing suspended 
graphene and graphene on BN EH EH in terms of line 
widths and of the variety of the observed excitations. 
The numerous features in Fig. |^b) dispersing with the 
magnetic field can be attributed to inter-Landau-level 
electronic excitations [13 m], and in the present exper¬ 
iment they can be observed up to energies as high as 
5000 cm“^. We observe B_n,n and „ (together 

with its electron-hole-symmetric counterpart L-n,n+i) 
with n up to 6 . Their evolution with magnetic field can 
be reproduced using the value v = ( 1.01 ± 0 . 01 ) x 10 ® m/s 
for the electronic velocity in graphene. 

The monotonous rise of the electronic excitation en¬ 
ergies with the magnetic field in Fig. [^b) is seen to be 
interrupted at several particular values of the magnetic 
field, identified by white and black boxes. These values 
of B match very well those corresponding to the reso¬ 
nances En = fy^KX for n = 1,2,3, shown in Table |T] 
In the classification of Sec. all these resonances corre¬ 
spond to case (hi) with ^ = 1 , so they are supposed to 
produce double anticrossings for L-n,n+i ex¬ 

citations, and triple anticrossings for L-n,n excitations. 
It is clearly seen in Fig. ib) that the resonance of L-n,n 
excitations with the two-phonon excitation occurs at an 
energy 2huJx ~ 2570 cm“^. This energy is different from 
that of the 2D band feature which involves phonons away 
from K point and which appears at higher energies (about 
2700 cm“^ for a 780 nm excitation). 

The difference between double anticrossings for anti¬ 
symmetric excitations and triple anticrossings for sym¬ 
metric excitations can be clearly seen in the Raman scat¬ 
tering spectra. In Fig. Qa),(b) we show the spectra cor¬ 
responding to the ± 0,1 and to the ±-i,i excitations, re¬ 
spectively, for selected values of the magnetic field close 
to B = 18 T. In Fig.|^a), the single peak observed at low 
magnetic field splits into two components with a grad¬ 
ual transfer of its intensity to the higher energy peak as 
El crosses the phonon energy. This effect on the ±o,i 
excitation has also been observed recently in infra-red 
magneto-transmission experiments performed on multi¬ 
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Figure 4 . a) Spectra of Lo,i electronic excitation in the 
range of magnetic field where it is in resonance with K point 
phonons, showing an avoided crossing behavior, b) Spectra of 
the ±-1.1 electronic excitation in the same range of magnetic 
fields as in a), and showing a triple avoided crossing. c,d) ±o,i 
and ±-1,1 energies, respectively, as a function of the magnetic 
field (black dots) together with calculated lines (solid lines). 


layer epitaxial graphene [28]. The evolution of the ±-i,i 
excitation in this range of magnetic fields is completely 
different, as seen in Fig. |^b): at the resonance, a three- 
peak structure appears. A similar triple avoided crossing 
can be observed also on the ±-2,2 when it is tuned to 
2hioji at B ~ 9 T, and when L-i 1 and L- 2,2 are tuned to 
2tiu>r, at B ^ 28 T and B ^ 14 T, respectively. Because 
the observed electronic excitations show a splitting, one 
can conclude that we are in the strong coupling regime 
and that magneto-polarons are formed at the resonant 
magnetic fields. 

More avoided crossing features can be seen in Fig. [^ 
where we present the evolution of the maxima of the ob¬ 
served electronic excitation Raman scattering features as 
a function of the magnetic field. In this figure, in addi¬ 
tion to the type (iii) resonances involving iF-point and 
F-point phonons (red and blue points, respectively) men¬ 
tioned in Table [Tj some resonances of the type (i) reso¬ 
nances involving K-point phonons can be seen (shown 
by orange points). They correspond to the condition 
En — El = huJK, and are expected to occur at B « 
12.3, 8.9, 6.9 and 5.6 T for n = 5, 6 , 7, 8 , respectively. 
We present in Fig. [^a) characteristic Raman scattering 
spectra between B = 9 T and B = 14 T and showing, 
among others, the ±- 4,5 excitation. Close to B = 12 T, 
the corresponding Raman feature splits into two compo¬ 
nents and then recovers its single-peak shape. Similar 
behavior is observed near B = 9 T on the ±-5,6 exci¬ 
tation. Our experimental resolution is not sufficient for 
a quantitative analysis of this data. We have not been 
able to resolve avoided crossings at the same values of B 






















7 



B(T) 


Figure 5. Magnetic field dependence of the main features ob¬ 
served in our Raman scattering spectra. Red points indicate 
resonances involving K-point phonons, blue point indicate res¬ 
onances involving F-point phonons and orange points indicate 
the effects of the Si _i) = Ek resonance. Green circles 
represent the case (iii) resonance with C — 2 observed on the 
L_ 3,4 excitation. The horizontal gray bar masks the residual 
contribution of the 2D band. 


on S_ 5,5 and S-6,6 features. A possible reason for this 
is that for L_n,n the splitting is necessarily accompa¬ 
nied by the intrinsic broadening, as discussed in Sec. EYl 
which makes it harder to resolve than for asymmetric 
„ excitations. No type (i) resonances with F- 
point phonons are observed either; this may be due to 
weaker coupling of electrons to these phonons than to 
K point phonons. Finally, the feature seen at B = 18.7 T 
on A_ 3,4 line corresponds to case (iii) with ( = 2. The 
Raman scattering spectra of this excitation are presented 
in Fig. 1^) for B = 15.5,17.5 and 18.5 T. The single 
peak observed below R = 15 T transforms into a broad 
feature at high magnetic fields. This evolution is also 
presented in Fig. in the form of large green circles. Be¬ 
cause of the lower sensitivity of our experimental setup 
at such high energies, we cannot fully resolve the line 
shape nor can we determine the number of components 
contributing to this feature at the resonance, although 
we expect that in this range of magnetic fields the reso¬ 
nance — El = El — Eq = hwK is relevant, so a triple 
avoided crossing should be observed. 

It is natural to describe a double or a triple avoided 
crossing phenomenologically by an effective two- or three- 



Raman shift (cm-1) 



Raman shift (cm'^) 


Figure 6. a) Raman scattering spectra (5 = 0 T sub¬ 
tracted) for selected values of the magnetic field showing 
the L_ 5 (_ 4 )_ 4 ( 5 ) excitation splitting around R = 12 T. At 
i? = 14 T, the principal electronic excitations are identified, 
b) Raman scattering spectra measured at R = 15.5,17.5 and 
18.5 T in the range of energy of the I/_ 3,4 electronic excita¬ 
tion. 


level problem with the Hamiltonian 


tj‘2 X 2 
-^eff — 


tt3x3 _ 


En + En±l C'iHv/Ib \ 

(3a) 

C'lhv/lB 

Endzl T EUph J 

2 En 

Cihv/lB 0 \ 


CiHv/Ib 

En+hujph C2hv/lB 

(3b) 

0 

C2fiv/lB SlkUph ) 



Here Wph is either ujk or wr, depending on which reso¬ 
nance is considered. The coupling between different lev¬ 
els is characterized by the natural energy scale hv / Ib and 
the dimensionless coefficients C), Ci, C 2 . For the mo¬ 
ment, we treat them as phenomenological fitting param¬ 
eters, different for each avoided crossing. Their values, 
extracted from the fit to the experimental data are given 
in Table [H] for the crossings which are sufficiently well 
resolved to allow for a quantitative analysis. 


The validity of the effective models (3a), (3b) is dis¬ 


cussed in detail in Sec. IV D where we also give approx- 



















imate relations between the phenomenological couplings 
C[,Ci,C 2 and the dimensionless electron-phonon cou¬ 
pling constants It is important that all reso¬ 

nances involving a given phonon (at either K oiT point) 
are described by a single electron-phonon coupling con¬ 
stant, so the analysis of the whole set of avoided crossings, 
besides producing an estimate for Aif, Ar, provides also a 
consistency check for the theory. The whole data is rea¬ 
sonably well described by \k = 0.05, Ar = 0.03, and the 
comparison between theoretical and experimental values 
of the coefficients C[,Ci, C 2 is given in Table |n] 

IV. THEORY 
A. The model 

In the Landau gauge, = —By, the electronic 
states in graphene are labeled by four quantum numbers: 

(i) Landau level index £, which is an integer running from 
—00 to -boo (in the previous sections, we used n = \£\)] 

(ii) the X component of the momentum p, which can be 
taken to run from 0 to Ly/l^ with spacing ‘1 'k if one 
considers a finite rectangular sample of the size x Lj,; 

(iii) the valley index, K or K'■, (iv) the spin projection, 
which does not play any role in the following, so it will be 
omitted. The energy depends only on the Landau level 
index, 

Ei = ^ ^/^\ sign^, (4) 

where the Dirac velocity v « 1.0 x 10® m/s « 6.6 eV • A, 
the magnetic length Ib « k)/yjB/{\ T), and we set 

h=l throughout this section. 

We consider electron coupling to two kinds of phonons: 
(i) the Aig transverse optical phonons near the K or 
K' points, whose states are labeled by the valley index 


K, K', and the wave vector q counted from the corre¬ 
sponding point (K or K')] (ii) the i? 2 g optical phonons 
near the T point, labeled by the wave vector q and one 
of the two circular polarizations tr = ±. We will neglect 
the phonon dispersion, so the frequencies of (i) and (ii) 
are wr « 196 meV and ujk ~ 160 meV. It should be 
emphasized that the most important contribution to the 
dispersion of phonons near the T or A points comes from 
their coupling to the Dirac electrons which produces the 
Kohn anomaly [Ml El] . The coupling to the Dirac elec¬ 
trons is included in the Hamiltonian below, so the Kohn 
anomaly should not be included in the phonon disper¬ 
sion to avoid double counting. Moreover, the spectrum 
of Dirac electrons is strongly modified by the magnetic 
field, so the resulting contribution to the phonon disper¬ 
sion is quite different from the one at H = 0. Thus, what 
is neglected here, is the mechanical part of the phonon 
dispersion which is due to the interaction with electrons 
in the a bands. 

The Hamiltonian is taken in the form 

H = Hll + Hph,K + Hph,r + (5a) 

The non-interacting part of the Hamiltonian is given by 

00 

Hhh = ^ ^ 

£— — oq P 

Hph,K = i (5c) 

q 

l?ph,r = (^q.+fcq.-e + > (5d) 

q 

where , c are the fermionic creation and annihilation for 
electrons in the corresponding single-particle states. The 
electron-phonon coupling is the same as in Ref. |M|) but 
written in the Landau level basis: 


Hk = 


T T 'A,f( q) e ® (^q,if4,p-g,„/2,K^^',P+9x/2.^' + ^q,if'4,p-g,„/2,K'^^'. P+9x/2.i<') +h-C., 

q 1,1',p 


Hr = 


\ T T q) e ^ (jyll,u^\,p-q„,/2,K^^',P+<I^/‘^,K ^q.o'^J,p-(jj,/2,if'd^tP+9x/2,if') +h.C.. 

V ci,a l,i',p 


(5e) 

(51) 


where “h.c.” stands for the Hermitian conjugate, and 
the dimensionless coupling constants Apf,Ar are defined 
as in Ref. [Mj- While the value Ar « 0.03 is established 
relatively well Uni 1321133] and is also in agreement with 
our experimental results discussed in the previous sec¬ 
tion, the constant Xk is subject to Coulomb enhance¬ 
ment [MIES], and thus is substrate-dependent. For our 
samples of graphene on graphite, the Coulomb renormal¬ 


ization effects have been shown to be strongly suppressed 
by substrate screening m, so the enhancement should 
be relatively weak; our experimental data are reasonably 
well described by Xk = 0.05 (see the previous section). 

The coefficients J£f/(q) = [,^|'/(—q)]* and (q) = 
q)]*, represent the matrix elements of the 
phonon-induced potential between the electronic eigen- 
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transition 

resonant levels 

C'l or Cl (exp) 

C[ or Cl (th) 

C 2 (exp) 

C 2 (th) 

Lo,i 

El = huJK 

0.044 ±0.004 

0.045 



-^-1.1 

El = huJK 

0.056 ± 0.003 

0.063 

0.084 ± 0.005 

0.070 

L-1,2 

3 

II 

Cb 

0.042 ± 0.003 

0.045 



T-1,2 

E 2 = hu)K 

0.049 ± 0.004 

0.045 



L-2.2 

E 2 = hoJK 

0.055 ± 0.004 

0.063 

0.088 ± 0.005 

0.072 

L-2,2 

E 2 = huY 

0.044 ±0.004 

0.049 

0.053 ±0.005 

0.057 

L-2,3 

E 2 = huJY 

0.046 ± 0.003 

0.035 




Table II. The dimensionless couplings C'i,Ci,C 2 in the effective two- and three-level models (Sal, ([Sb]) for various avoided 
crossings observed in the experiment. The theoretical values are obtained from relations given in Sec. |IV D| with dimensionless 
electron-phonon coupling constants Xk = 0.05, Ar = 0.03. 


states in the magnetic field, and are given by 

JIM) = V(i + M(i + <5t'o) - 

- (1 - ^m)(l - h’o) ^ 

(6a) 

Jtti^) = (1 - <Jm)sign£ ^ ° J|^|_i,|^/|(q), (6b) 

'4>(q) = (1 - <J£'o)sign/ ^ J\i\,\i'\-M)- (6c) 

The formfactors, 


/ \ _ ( _1 \m-|-min{n,m} I Cnin{Tl, ? 7 l}! 

> \/max{n,m}! 


X e 


- 9 ^^s /4 f + iQy 


\n—m\ 


qlB\ 

9 / xV^J 

xLiirt:.'it(9^^i/2)=^;;,n(-q), (7) 


min{n,m} 

are defined via the associated Laguerre polynomials, 


k=0 


k\ {m — k)\ + k)\ 


( 8 ) 


which satisfy the following orthogonality relation: 


OO 


the same way. The spectral function can be obtained as 
— (I/tt) Imn„ j (r2) from the excitation propagator 


n„.j(fl) = -i 


2 Trlg 

LxLy 


T7?n,,W<,(0))e 


Aflt 


dt. 


( 11 ) 

where the average is over the ground state and T stands 
for the chronological time ordering. 

The excitation propagator (111 is calculated using the 
standard zero-temperature diagrammatic technique [36], 
which is constructed from the following basic elements. 


• The electron Green’s function, represented by a 
solid line, carries the energy variable e and the in¬ 
dices of the single-particle states. It depends only 
on the energy and the Landau level index: 


GM 


1 

e — E(^+ f 0+ sign £' 


(12a) 


The sign of the infinitesimal imaginary part fO’*' 
corresponds to the Landau level filling: the lev¬ 
els with £ < 0 are filled, those with £ > 0 are 
empty. On the level £ = 0 only a fraction /g of 
states is assumed to filled. This can be modeled by 
introducing a fictitious dependence of the infinitesi¬ 
mal imaginary part on the momentum p: we choose 
at random a fraction /g of the momenta for which 
we set Gp(e) = l/(e — iO’*'), and for the rest it is 
Gp(e) = i/(e + i0'^). In the subsequent calculations 
it is equivalent to simply setting 


We are not interested in the Raman process itself, 
which was studied in detail in Ref. [22]. Our main fo¬ 
cus is the spectral function of the electronic excitation, 
which was created in this process. The creation operators 
for the electronic excitations observed in the experiment, 
those of the type —n —>■ n, — (n ± 1) —n, can be repre¬ 
sented as 

~ ^n.v.K^-{n+j),p,Kj j = 0,±1, (10) 

P 


Gfco(e) 


l-/o 
e -I- *0+ 


fo 

e — 10+ ’ 


(12b) 


• The phonon Green’s function, represented by a 
wavy line, carries the frequency variable uj, the 
two-dimensional momentum q, and one of the four 
labels K, K\+,— (the latter two correspond to 
the circular polarizations of the phonons near the 
r point). The Green’s function is given by 


where we have arbitrarily chosen the states in one valley, 
K, as the states in the other valley behave exactly in 


7?if,r(w) 


1 

u! — ujK,r + *0+ 


(13) 
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(a) 



(b) 


1 . Asymmetric transitions 

For the asymmetric transition — (n ± 1) —>■ n the res¬ 
onant contribution to n„_±i(r2) is given by the sequence 
of diagrams shown in Fig. I^a). It is sufficient to dress 
the electronic Green’s function G„(e) by the self-energy 


o+dD-^c: 


+ ... 





E^{e) = \KV^J 


cPq 




(27r)2 e — luk — En' + i0+ 


i/2 


S„(e) = Aru 


e — lok ~ En' + *0+ 

l4+„,(q)P + |J-„,(q)P 


(27r)2 e — wr — En' + t0+ 




T/2 _ 


6 — (x^p — Eji! -\- 

ApcrU^ 


47rl| 


so that the dressed Green’s function is given by 


(15a) 


(15b) 

(15c) 


Figure 7. Diagrams contributing to the excitation propagator 
in the case (i) of Sec. (a) for the asymmetric transition 
L-(n±i),n, (b) for the symmetric transition L-n,n, with the 
vertex function given by (c). The thick solid line represents 
the dressed electronic Green’s function G±n{e), the shaded 
box represents the vertex function. 


(e) - S±„(e) 

^ _ e T (^ic,r + En') _ 

[e=F (^n - jO+)][e=F (^« - *0+)]’ 

where E^ are the two roots of the equation 


We do not include the negative-frequency part of 
the phonon Green’s function as the resonant ap¬ 
proximation will be used in all subsequent calcula¬ 
tions. 


• The electron-phonon vertices are given by the co¬ 
efficients in the interaction Hamiltonians Hk.v in 
Eqs. (l5e|), ISfl). 


In the absence of interactions, the excitation propaga¬ 
tor is given by 


(E-F;„)(E-<-ccip,r)-F|,r=0, (17) 

representing the energies of the two hybrid (magnetopo- 
laron) states. The resulting excitation propagator, 

n„,±i(H)= y|iG„(e)G_(„±i)(e-H) = 

_ Eji^i E-fi' Q\ 

“ (O - En±i - E+){n - E„±1 - EnY 

where the positive infinitesimal imaginary part of H is 
omitted for compactness, coincides with the upper left 
element of the matrix (H — E[^g^)~^, where 


nS(f^) = / G„(e) G_(„+,)(e - H) = 


~ (En -|- En+j) -|- *0+ 

SO the spectral function is given by just a single (5-peak 
6 {n-{Er,+Er,+,)). 


H. 


eff 


En±l + En 
Vk,T 


En±l 


Vk,t 

+ En' + UlKX 


(19) 


can be viewed as the Hamiltonian of an effective two-level 
model. The resulting spectral function is given by the 
sum of two peaks 5{^l — En) the weights determined 

by the corresponding eigenvectors. 


B. Irrational case 

Let one of the phonon frequencies ujk,t be close to 
the difference En — En' and C is irrational (equivalently, 
y/njvf is irrational). This is case (i) of Sec. |h[ 


2. Symmetric transitions 

As discussed in Sec.|TTj in the case of a symmetric tran¬ 
sition —n —>■ n, it is natural to expect the problem to 
be analogous to that of three coupled levels. However, 
the situation turns out to be qualitatively different for 
phonons near K,K' points and for those near F point. 
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For the resonance En — En' ~ ojk, it is sufficient to 
use the dressed Green’s function both for the electron on 
the level n and for the hole on the level —n, Eq. (16), 
and to calculate as in Eq. (14), but using the 

dressed Green’s functions instead of the bare ones. This 
corresponds to the first term in the sequence of diagrams 
shown in Fig. [^b), and gives 


i^-K- E+m -K- E-) - Vl 
{Sl-2Et){Sl-2En){Sl-Et-E^)' ^ 


which is indeed equivalent [in the same sense as above, 
i. e. matching the upper left matrix element of (fl — 
ideff)”^] to an effective three-level system with the Hamil¬ 
tonian 


E. 


2E, 

3 ^ 3 = ( V2Vi 


eff 


0 


0 

K En-\- En' -f LOk Vk 


V2Vi 


K 


( 21 ) 

Accordingly, the spectrum consists of three 5 peaks. 

For the resonance En — En' ~ WP) there are vertex 
corrections. This gives rise to the rest of the terms in 
Fig. [3b) , whose contribution is of the same order as the 
first one. The vertex, in turn, is represented by an in¬ 
finite sum of diagrams, whose first terms are shown in 
Fig. E c). Physically, it corresponds to the electron and 
the hole exchanging the emitted phonons before their re¬ 
absorption. This process is possible only for the F point 
phonons; for the K, K' phonons, those emitted by the 
electron and by the whole necessarily belong to differ¬ 
ent valleys, and cannot be exchanged, so all diagrams in 
Fig. [3b) vanish except for the first one. 

Inclusion of vertex corrections significantly increases 
the level of computational difficulty. Even though the 
diagrams for the vertex function can be summed, this 
summation reduces the problem to a system of inte¬ 
gral equations for functions of momentum q with a non- 
separable kernel. These equations are given explicitly 
in Appendix [Aj where a different but equivalent formu¬ 
lation of the problem in terms of magnetopolaron wave 
functions is given. Solution of the resulting equations is 
beyond the scope of the present paper, but the qualita¬ 
tive properties of the result are worth discussing. 

The dramatic effect of the phonon exchange is that the 
problem can no longer be effectively described in terms of 
a few coupled discrete levels. Let us note first that such 
description arose in the previous case only because we 
neglected the phonon dispersion as well as the dispersion 
of electronic excitations as a function of the total momen¬ 
tum. Had we included the dis persi o ns, al ready in the sim¬ 
plest case of the self-energies (15a), (15b I, they would not 
reduce to a simple single-pole form, so the spectral func¬ 
tion Imn„_±i(H), instead of being a sum of two 5 peaks, 
would broaden into a continuum formed by excitations of 
different momenta. In the case of phonon exchange, even 


if one originally starts from non-dispersive excitations, 
the dispersion is produced as a result of the momentum- 
dependent coupling between the one-phonon and the two- 
phonon sectors, so the spectral function also broadens 
into a continuum with a non-trivial shape. It is impor¬ 
tant that this broadening should be of the same order as 
the splitting at resonance, as both are determined by the 
same energy scale Hr- 


3. Phonon dispersion 


The phonon frequencies acquire a resonant correction 
when ujx or wp is close to the energy i?„ — E_n' of a 
transition between filled and empty states. For ^njn' 
irrational, such resonances are distinct from those anal¬ 
ysed above. The shifted frequencies can be found from 
the equation w) = 0, where the dressed phonon 

Green’s functions ^re given by the diagrams 

shown in Fig. [3a). Summation of the series gives 


(q, w) = w - Wif - 4V^ ■ 


uj — Efi — Efi' + zO"*” 


(22a) 


n-l. X + l'^n',-n(q)P 

Up (q, uj) = to — LOr — oVp 


CU — En — En' 


(22b) 


The coefficient is twice larger for the F phonons (8 in¬ 
stead of 4) because they can excite intravalley electron- 
hole pairs in any valley, while a K phonon can excite an 
electron in K' valley and a hole in K valley only (the 
spin degeneracy has also been taken into account). 

The roots of the quadratic equation Z))^^p(q, w) = 0 
for each q give two frequencies of the coupled excita¬ 
tions which are superpositions of the phonon and the 
electronic excitations. As J„„'(q) quickly falls off for 
q 3> 1/^B, the resonance modifies the phonon dispersion 
in the region q ~ 1/?b, \En + En — u]K,T\ ^k.t, and the 
typical magnitude of the frequency shift in this region is 
given by Vkx ~ \/^kx v/Ib- This square-root depen¬ 
dence on the coupling constant should be contrasted to 
the Kohn-anomaly correction to the phonon dispersion 
without magnetic field, which is of the order of Xkx'^ 9 - 


C. Rational case 

The smallest Landau level indices n, n' which produce 
a non-integer rational C are n = 9, n' = I, leading to (^ = 
3/2. Even this resonance, Eq — Ei = ujkx^ beyond our 
experimental resolution. Still, it is not too far from the 
observed resonance Eq — Ei = ujkx^ so the perspective 
of observing the Eg — Ei = ujkx resonance in the future 
is not totally hopeless. For more complicated cases, such 
as Eiq - El = UJKX giving C = 4/3, or E 25 - Eg = ujrx 
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(a) 



(b) 


+ 


<£3 


(c) 


o = CS + 


+ 


Figure 8. Diagrams contributing to the excitation propaga¬ 
tor n9,±i(f2). (a) Dressed phonon Green’s function, (b) Ex¬ 
citation propagator for the asymmetric transitions L_g,g, 
L- 10,9 with the electronic self-energy given by (c). The thick 
wavy line represents the dressed phonon Green’s function 
w), the shaded circle represents the electronic self¬ 
energy. 


giving C = 5/2, there is not even such hope. Thus, here 
we focus on the simplest case Eg — Ei = ujkx- 

The difference of this resonance from those discussed 


in Sec. IV B is that the condition Eg — Ei = cok.f au- 
tomaticaily yields ujkx = Ei — E-i. Thus, the emitted 
phonon can be resonantly reabsorbed by an electron on 
the level n = — 1 , creating a second electron-hole pair, in 
addition to the one initially produced by the photon (the 
energy of the resulting excitation is the same as that of 
the original excitation, as Eg + E-i = 2Ei). Thus, we 
dress the phonon Green’s function, as shown in Fig.|^a): 


3(q,w) = 

i?r(q,w) = 


CJ — (jJk — 


UJ — (jJY 


UJ - 2Ei + i0+ 

8V/g|Jf.-i(q)P 

UJ - 2Ei + f0+ 


-1 


-1 


(23a) 

(23b) 


For the asymmetric transitions, T-s.g and L_io, 9 , the 
excitation propagator ng_±i(r 2 ) is given by the series in 
Fig. ^b), where the electronic Green’s function is dressed 
by the self-energy insertions. As in Sec. |IV B 2\ there is a 
difference between the coupling to AT, K' phonons and to 
F phonons, related to exchange. Namely, when the initial 
excitation (say, T- 8 , 9 ) is resonantly converted into T_ 8 ,i 
and one K or K' phonon, and then into L_ 8 ,i+T_i 1 , the 
two electrons on the n = 1 level must belong to different 
valleys, and thus are effectively distinguishable. As a 
result, in the series for the self-energy, shown in Fig.j^c), 
only the first term survives for K, K' phonons, resulting 


in a relatively compact expression: 

Sf(e) = Xkv^J ^ |J9\(q)3K(q,e- Ai) = 

00 

^ [ V^{e-Eg){^-6)^ee-^d^/{2-9\) 

J {e-Eg){e-uJK-Ei)-V^{^-2)^e-i- 

0 

(24) 

Then the excitation propagator is straightforwardly eval¬ 
uated as 


n9,±i(fi) = 


1 


n-{Eg+ Eg±i) - Sf (fl - A 9 ±i) ■ 


(25) 


As discussed in Sec. |IV B 2[ as soon as excitations in¬ 
volved in the decay of the initial electron-hole pair ac¬ 
quire a dispersion, the peaks in the spectral function 
broaden and the situation can no longer be reduced to a 
simple effective picture of a few coupled discrete levels. 
This is what happens here, due to the phonon dispersion 
acquired from the resonant coupling to A- 1,1 electronic 
excitation, as seen from Eqs. (23a), (23bI. However, in 


the present case, this effect turns out to be strongly sup¬ 
pressed due to the large Landau level index n = 9. Let 


us formally expand the self-energy Eg (e) from Eq. (24) 
at V^/[{e - Eg){e -ujr - Ex)] 00 : 




(e) = 


1/2 


e-ujK - El 

5 


1 + 


1/2 


2048 (e — £'9)(e — ujk — 


(26) 


At resonance, ^ {e — Eg){e — ujk — Ei), the second 
term in the square brackets is still small because of the 
numerical factor. The origin of this numerical smallness 
is that the numerator of the integrand in Eq. (24) has a 
maximum at ^ « 12 and is very small at ^ < 7, while the 
function (^ — 2)^e“^ in the denominator is quite small 
for ^ > 2. Physically, this means that the electronic 
transition from level n = 9 to n = 1 is accompanied 
by emission of phonons with wave vectors q = \/2^jlB 
which are too large and do not efficiently couple to A- 1,1 
excitation. 

In the first approximation, if we neglect all terms in 
Eq. (26) except the first one, which is equivalent to re¬ 


placing the dressed Green’s function A>if(q, oj) by the 
bare one, i. e., neglecting the coupling of the phonon to 
A_i^i excitation, we naturally obtain the effective 2x2 
problem. 


tj2X2 


Eg Vk 

Vk El + UJK 


(27) 


whose two eigenvalues Eg determine the positions of two 
6 peaks in the spectral function. When the full self¬ 


energy (24) is taken into account, the two peaks (i) shift 
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Figure 9. The dimensionless spectral density, 

— (Fk/tt) Imn 9 ,±i(fl), versus the dimensionless frequency 
2 = {yi — Eq/2 ^ E-q±-i) /Vk for the case of exact reso¬ 
nance, Eg — El = ujk- The square-root-type singularities 
at z = ±2e“^ ~ ±0.27 correspond to the maximum of 
— 2)e“^^^, the dimensionless energy splitting between the 
phonon and the excitation [see Eq. ( |24[ )]. The spectral 

density for \z\ > 2e“^ is extremely small, ~ 10“®, but finite; 
it vanishes exactly for \z\ > 2. The two main peaks are 
located near z = ±1 (not shown). The singularity at 2 —>■ 0 
is logarithmic and comes from large 


(the shift being equai in magnitude and opposite in sign 
for the two peaks), (ii) broaden [as ImE|^(i?^) ^ 0 , ai- 
though very smaii], and (iii) iose a smaii part of their 
spectrai weight. This part of the spectrai weight is trans¬ 
ferred to the continuum concentrated in the vicinity of 
the energies of the uncoupied excitations; in this region 
the spectrai function has a pecuiiar shape, as iiiustrated 
in Fig. I^for the case of exact resonance. Eg — Ei = ojk- 
At the same time, the two main peaks have neariy 
Lorentzian shapes, determined by the poies of Ilg j(f2) 
in the compiex piane. The poies are iocated at 


— Eg^j ± 


Eg + El + UJk 
2 


± ( 28 ) 

where u± « 1.00251 ± 0.66 x 10“®i are the solutions of 
the equation 


0 = E{u) = 1 - 


/ 


(g-6)^g^e ^ 

u — (^ — 2)2e“5 


JL 

2-9!’ 


(29) 


and the sum of the residues at the two poles is given by 


dE{u) 

du 


0.9975. 


Thus, corrections to the simple two-level picture of 
Eq. (27) are quite small and will be masked by other 
effects, not included in the present model. 


The same can be said about the case of coupling to 
r point phonons. Even though the self-energy, corre¬ 
sponding to the sum of diagrams in Fig. 


±(«) = E E 


(i2qi...d2q^ ~ 

\k *^91 (^l) ['^91 i^k)] X 


(V^l) 


k—1 (Ti ...(Tfc 

k-1 

1=1 


X Vt 


2 k - Pic(qi; £ ~ El) . . . DKjfjkii- El) 

{Eg - e - 


(30) 


is not easily evaluated, the factor (q) will again result 
in an efficient coupling to phonons with relatively large 
wave wave vectors, which disperse very little. The same 
argument applies to the case of the symmetric transition 
L-g^g. Thus, for all practical purposes, the case of ra¬ 
tional (C can be considered equivalent to the irrational 
case. 


D. Integer case 


Among integer values of C, the most interesting one 
is = 1 , that is, the resonance between the excita¬ 
tion En — Eq and a phonon. This corresponds to the 
most pronounced features in the experimental data pre¬ 
sented in Sec. |III| The next value, C = 2, corresponds 
to the resonance — En /4 = u}K,r ■ A signature of this 
resonance is seen as broadening of the ±- 3,4 transition 
around S « 18 T, but the resolution is not sufficient for 
a reliable quantitative analysis of the data. At the same 
time, the experimental data for C = 1 resonances do allow 
for a quantitative analysis. 

From the theoretical point of view, description of reso¬ 
nances with integer represents an extremely hard prob¬ 


lem. In Fig. 10 we show diagrams contributing to the 
excitation propagator n„^±i(n) in the simplest case of 
the asymmetric transition T_(n±i),n when — Eg is 
resonant with a phonon (C = 1 ) and n > 1 [for Lq.i 
transition, the diagrams for are the same as in 

Fig. [Io][c)]. The role of electrons on the partially filled 
n = 0 level is emphasized by representing their Green’s 
function by a dashed line. In fact, we are unable to 
identify a dominant sequence of diagrams. The physi¬ 
cal reason for this is the macroscopic degeneracy of the 
electronic non-interacting ground state due to partial fill¬ 
ing of the n = 0 Landau level. As a result, conversion 
of the initial electronic excitation with energy i?„ into 
the phonon and back can be accompanied by emission 
of an arbitrary number of intra-Landau-level excitations 
which cost no energy and can be exchanged in arbitrary 
order. Thus, various diagrams can be generated. In the 
resonant region, \En — ojk,t\ ~ Vk,t, ah these diagrams 
are of the same order, and there is no parameter for se¬ 
lection of a treatable subset of diagrams. As a result, the 




















14 


(a) 



where 0 < /o < 1 is the filling of the zero Landau 
level. For the resonance with F point phonons, the 
phonons corresponding to the two circular polarizations 
are dressed differently (unless /o = 1/2), so the above 
equations should be modified by (i) replacing Dk —>■ 

J* —>■ J'^, and (ii) introducing an additional factor of 
two in the first equation, —>■ 8V^, due to the valley 

degeneracy. Restricting ourselves to the electron-hole- 
symmetric case with /o = 1/2, we eliminate the phonon 
Green’s function, and obtain a single algebraic equation 
for G„{e) = -G_„(-e): 


(d) 


0 




-I- 


-I- 



Figure 10. (a) Diagrams contributing to the excitation prop¬ 
agator n„,±i(D) for n > 1. (b) The first terms of the series 
for the electronic self-energy which defines the dressed Green’s 
functions G±„(e). (c) Dressed phonon Green’s function, de¬ 
fined via the phonon polarization operator, the first diagrams 
for which are shown in (d). Thick solid lines represent the 
dressed electron Green’s functions on the nth Landau level. 
Dashed lines represent the bare electron Green’s functions on 
the partially filled zero Landau level. Thick wavy lines rep¬ 
resent the dressed phonon Green’s function. Shaded circles 
represent the electronic self-energy. Shaded rhombi represent 
the phonon polarization operator. 


coupled levels broaden into a continuum, and we are not 
even able to write a closed set of integral equations to 
describe this broadening. 

A qualitative description of this broadening effect can 
be done by keeping only the first diagram in Fig. @(b) 
and the first two diagrams in Fig. @(d). This proce¬ 
dure, analogous to the GW approximation for interacting 
electrons (see Ref. m for a review), gives a closed system 
of self-consistent equations for the dressed electron and 
phonon Green’s functions, G and D. For the coupling to 
the K point phonons this system has the following form: 


b-^\ci,uj)=uj-uJK-W^ [/o|J^o(q)PG„(a;) 


- (l-/o)|JoVn(q)l'G_„(-a;)' , 

(31a) 

G'„^(e) = e- En - {1- /o)Aiyu^x 


(31b) 

G_n(—e) — —e + En /oAku^x 


(31c) 


G„(e) 


= e-A„- 


<Vl|./no(q)Pd^q/(27r)^ 
i-ujK- 2y|| J„o(q)pG„(e) 


• (32) 


The equation for the resonance with the F phonon is 
obtained from this one by replacements lok —t wp, 
l-^«o(q)P |-^n-i,o(q)P- It is convenient 

to introduce dimensionless variables z = (e — ujk)IVk 
and g(z) = VxGnie), as well as the dimensionless de¬ 
tuning S = {En — (jJk)IVk- Then the Raman spectrum 
for the asymmetric excitation, (—I/tt) Imn„_±i(n) = 
(—I/tt) ImG„(r2 — i?„±i) oc (—I/tt) Im 5 (z). The lat¬ 
ter quantity, obtained from the numerical solution of 
Eq. (32) (see Appendix for discussion of some prop¬ 
erties of this equation) is plotted in Fig. 0 for n = 2 
and two values of detuning <5 = 0,1. 

The same discussion can be applied to the magne¬ 
tophonon resonance in the Raman scattering on the F- 
point phonons, seen as the resonant modulation of the 
G peak mm- The theory developed in Refs. uma, 
when considered close to each resonance, essentially maps 
on an effective two-level problem. This is correct for 
transitions not involving the zero Landau level, whose 
partial filling makes the picture for the Lop transi¬ 
tion more complicated, on the same level of difficulty 
as discussed above. Indeed, the Raman spectrum in 
this case is proportional to the phonon spectral func¬ 
tion, (—I/tt) ImL)ip(q = 0,a;), given by the diagrams of 
Fig.[^b). The effect is illustrated in the same GW ap¬ 
proximation in Fig. Mh), where we plot the spectrum 
{-l/-K)lTa[l/{z - g{z))\ cx (-I/tt) Imi)r(q = 0,a;), ob¬ 
tained from the numerical solution of Eq. (32), modi¬ 


fied to describe the coupling to the F phonons, using 
hr as the unit of energy and introducing the shifted fre¬ 
quency z = [lo — ujy)IVt, and for two values of detuning 
(5=(L;i-wr)/Vr = 0and 1. 

The GW calculation, presented above, is not justified 
by any small parameter; near the resonance, the diagrams 
in Fig. [To](c,d) which are neglected, are of the same or¬ 
der as those that are retained. So there is no reason 
to believe that the complicated spectral shapes shown 
in Fig. 11 have any quantitaive significance. (For the 
same reason we do not study the symmetric excitations, 
L-n,n in the GW approximation.) Fig. 11 serves only 


as an illustration of the key fact: even in the simplest 
model used here, which neglects many effects (such as the 
dispersion of phonons and inter-Landau-level electronic 
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Figure 11. The dimensionless spectral densities (a) of the 
asymmetric electronic excitation, —(Vk/tt) Imn 2 ,±i(ll), and 
(b) of the r point phonon, (—Vt/tt) Im Ilr(q = 0,a;), versus 
the dimensionless frequency (a) z = {Q — lok + E-^±i)/V k, 
and (b) z = (uj — uJr)/Vr, as obtained from the GW approxi¬ 
mation. On both panels, the solid line correponds to the case 
of exact resonance, 5 = (l? 2 ,i — ojk.t) /VK,r = 0, the dashed 
line to 5 = 1. The spikes for <5 = 1 do not correspond to 
any divergencies in the spectral function, they are sharp but 
finite. 


excitations, their broadening by disorder, fine structure 
of intra-level electronic excitations on the zero Landau 
level), the electron-phonon interaction involving the par¬ 
tially filled Landau level broadens the discrete peaks into 
non-trivial spectral shapes. The details of these shapes 
are likely to be masked by the above-mentioned effects 
anyway. 


Still, we need a practical recipe for a quantitative anal¬ 
ysis of the experimental data. In Sec. Ill we have seen 
that the experimentally measured magnetic field depen¬ 
dence of the positions of the Raman peaks’ maxima is 
reasonably well described by an effective two- or three- 
level model, Eqs. (3a), (3b). To fix the effective param¬ 
eters Cj, Cl, (72, we match the level shifts obtained from 


the effective models (3bI, (3a) away from the resonance, 
\En — uJK,r\ ^ to the expressions for the energy 


shifts obtained perturbatively from the original electron- 
phonon Hamiltonian of Sec. |IV A 

First, let us find the correction to the energy 2En of 
the initially excited electron-hole pair, 

\sh} = cl^pj^C-n,p,K\0}, (33) 


where |0) is one of the many degenerate ground states 
with partially filled zero level. It is specified by occu¬ 
pations (0 or 1) of states with different momenta, the 
fraction of filled states being /q. Strictly speaking, the 
optically excited state is given by the symmetric linear 
combination of states (33) with all momenta p and in 
both valleys. However, the energy shift depends neither 
on p, nor on the valley, so we do not write the whole linear 
combination for the sake of compactness. The second- 
order perturbation theory in electron-phonon coupling 
gives 


SEeh = 


|Jno(q)P ^ 

(27r)2 E„ - UK E„-ujk’ 


(34) 


for the coupling to the K point phonons. For the T point 
phonons one should replace Vk Vt, uk —>■ wr- At 
the same time, the shift from the effective model (3b) 
is {Civ/IbY/{E n-UKp), which fixes (7i = i/A/fy/(47r) 
for symmetric transitions. The shift (34) for the symmet¬ 
ric excitation comes from the coupling to states with one 
phonon, which can be emitted either by the electron or 
by the hole. For asymmetric transitions only one of these 
processes is allowed. Thus, the energy shift is given by 
the same Fq. ( [M| ), multiplied by 1 — /o if the initial exci¬ 
tation involved level n, and by /o if it was —n, so the cou¬ 
pling matrix element of the effective two-level model (3a) 
is C[ = y(IA77())XRrj7(^ or C[ = h\KX/ 


To find the effective matrix element V 2 , we have to find 
the energy shift in the one-phonon sector. The states in 
the one-phonon sector are characterized by the phonon 
wave vector q and have the form 


|e,q) = (35a) 

|e,q,CT) = (35b) 


for the K and T point phonons, respectively. The states 
|0e) represents a ground state with partially filled zero 
level, but it may be different from the state |0) be¬ 
fore the initial excitation. In addition to states (35aI, 
(35b), which correspond to phonon emission accompa¬ 
nied by the hole annihilation, there are also states cor¬ 
responding to electron annihilation obtained by replace¬ 


ment c, 


n,p,K 


C-n,p,K, and of the label “e” by “/i”. 


The second-order correction to the the energies of 
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states (35a I, (35b I is given by 


T/2 


[(l-/o)-|Jno(q)|"], 


(36a) 


SE, 




e,q,+ 


En — UJr 


[(l-/o)-/o|Jn-i,o(q)P] , (36b) 




SEe,^.- = [(1 - /o) - (1 - /o)| Jn-i,o(q)n . 

(36c) 


The first term in the square brackets comes from coupling 
to the two-phonon sector by annihilation of the remaining 
electron, which requires an empty state on the zero level. 
This gives the factor 1 — /o, which should be replaced by 
/o for the corrections SEh^q, SEh^q^a-- The second term 
is due to coupling to the zero-phonon sector, when the 
phonon can be absorbed either by producing a hole on 
the level —n or a second electron on the level n. For the 
K point phonon both these process are allowed, so the 
occupation /g cancels out, while for the T point phonons 
either one or the other process is allowed for each of the 
two circular polarizations cr = ±. Obviously, for the 
corrections SE^.q, SEh.q.g the second factor is the same 


as in Eqs. (36a)-(36c). 


The energy shifts are different for states with different 
phonon wave vectors q, which is another manifestation 
of the effect that has already been discussed: the energy 
spectrum of the coupled system does not reduce to a few 
discrete levels, but is continuous. Thus, we have to decide 
which energy in this continuum should be associated with 
the discrete state of the effective three-level model (3b I. 
Clearly, some values of q are more representative than 
others, as they have larger weight in the Raman spec¬ 
trum. This weight is determined by the overlap of the 
exact eigenstates with the initial \eh) state. Away from 
the resonance, it can be determined in the first order of 
the perturbation theory in the electron-phonon coupling. 
For the coupling to the K point phonons, the weight of 
the states |e,q), |/i, q) is given by 


|(e/i|e,q)p = /o 


27rZ| Fl|Jno(q)|" 


LxLy {En — IOkY ’ 


|(e/i|/i,q)p = (1 - /o) 


27rZ| V^|J„o(q)|' 


ExLy {En a^ic)^ 
while for the F point phonons we have 

2Trl% F[?|J„_i,o(q)P 


|(e/i|e,q,-)| = fo 


ExEy {En ) 


ExEy {En cup)^ 


\{eh\h,q,-)f = \{eh\e,q,+)f = 0. 


(37a) 

(37b) 


(37c) 


1 / u\u I \|2 /I f \ V[?| J„_i o(q)P /'qvjn 

\{eh\h,q,+)\^ = {1- (37d) 


(37e) 


We choose the “center of mass” as the representative en- 


transition 

resonance 

Cl 

C 2 

L-i^i 

El = UK 

0.282 VAif 

0.315 VAk 

L-i^i 

El = LJr 

0.282 V/w 

0.315 VT/ 

L- 2,2 

E 2 = UK 

0.282 VAif 

0.323 VAk 

L- 2,2 

E 2 = wr 

0.282 VA/ 

0.331 VA/ 

L-3,3 

E 3 = Uk 

0.282 VAjf 

0.327 VAk 

L-3,3 

E 3 = ojr 

0.282 ^/A/ 

0.335 VA/ 


Table III. The dimensionless couplings Ci, C 2 in the effective 
three-level model (3b I for the three lowest symmetric elec¬ 
tronic transitions T- 1 , 1 , T- 2 , 2 , 7/-3,3 and K and F point 
phonons at half-filling of the zero Lan dau level (/o = 1/2). 
For the effective two-level model (Sal, we then have C( = 
Ci/V^. 


ergy of the continuum: 


^ llq^Ee,q\{eh\e,c^\'^ + Y.q^Eh.q\{eh\h,q)\'^ 
Otj = — 


En — 


5E+ = 




EJ(e%,q)P + EJ(e/i|/i,q)P 

(2n)! 

22"+i(n!)2j ^ 

(2n- 2)! 


5E_ = 


En — Wp 

(1 - h)V^ 

En — Wp 


2 / 0 ( 1 -/o)- 
1 - 


22 n-l((„_ l)!)2j ’ 


1 - 


(2n- 2)! 


22’^-i((n- 1)!)2_ 


(38a) 

(38b) 

(38c) 


These shifts should be m atched with that obtained from 
the effective model (3bl, {C 2 — Cl){v/I bY/{E n — 


The resulting values of C 2 for the first three transitions 
at half-filling (/g = 1/2) are given in Table m They are 
compared to the experimental values in Table |n] (end of 

Sec.|!TT|. 


V. CONCLUSIONS 

To conclude, we have presented a complete theoretical 
and experimental picture of magneto-phonon resonances 
in graphene, as observed in the magneto-Raman scat¬ 
tering spectra. By considering the case of an electronic 
inter-Landau-level excitation decaying into a phonon and 
another electronic excitation, we have extended the com¬ 
monly accepted picture of magneto-phonon resonance in 
graphene (i) to both F and K point phonons and (ii) to 
multiphonon processes. Such multi-excitation final states 
have numerous possibilities to fulfill momentum conser¬ 
vation, including intervalley electronic excitations. We 
have derived a simple classification of such magneto¬ 
phonon resonances and we have described the properties 
of each different type of resonances. In particular, res¬ 
onant splitting of three crossing excitation branches is 
observed experimentally for the first time. Such multi¬ 
phonon processes are particularly important for the un¬ 
derstanding of hot carrier relaxation in graphene. 

We analyze the experimental data phenomenologically 
in terms of effective two- or three-level models, and then 
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give a detailed discussion of validity of such effective de¬ 
scription starting from a microscopic theory in the reso¬ 
nant approximation. We also highlight the richness of the 
physics associated with the partially filled zero Landau 
level and its low-energy electronic excitations. When this 
level is involved in a resonance, the possibility to excite 
an arbitrary number of intra-level excitations with zero 
energy (i) leads to a genuine many-body problem which 
is quite hard to analyze theoretically, and (ii) results in 
complicated spectral shapes of the split peaks. 
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Appendix A: Magnetopolaron wave functions 

Let us first consider the case of coupling to K, K' point 
phonons, which results in Eq. (20). Let us look for the 
eigenstates in the form 

Id/) = d/^^le/r) + ^ d/^|e, q) + ^ q) + 

q q 

+ ^qi.q2lqi-q2), (Al) 

qi.q2 


the overall valley symmetry. The optical excitation cor¬ 
responds to the symmetric combination of states from 
the two sectors. 

The Schrodinger equation — flld/) results in the 

following system of linear equations for the coefficients in 
Eq. ([Ml: 

Xkv^ 

= -J2 + E (A3a) 

q q 

{n-En- - lok) = 

(A3b) 

(A3c) 

yS(L!-2E„,-2a;K)d'q,,q. = 

- (A3d) 


where the basis states are defined as 


Let us look for the solution in the separable form 


= \ E cl.p.KC-n,p,K|0), (A2a) 

V p 

V p 

X SE|0), (A2b) 

V ^ ^ p 

xbl^,\0), (A2c) 

V ^x-^y p 

X ^qi ,kKi 2 .k' 10), Q = qi -f q 2 • (A2d) 


Here |0) is the ground state. Note that, generally speak¬ 
ing, d/q^^qj ^ d'q 2 ,qj as the phonons belong to different 
valleys and thus are effectively distinguishable in the two- 

phonon state |qi, q 2 ). _ 

Eigenstates of the form ( |Al[ ) belong to the sector where 
the initial electron-hole pair \eh) is in the K valley. There 
are also eigenstates in the complementary sector, ob¬ 
tained by exchanging K -H- K' in the above expressions, 
with exactly the same energies, as the Hamiltonian has 


'I^q = -1 


I XkV^ 


LxLy 


Vk 




xT/'/ = 
q 


Kn'(q)]* ,j,h 
Vk 


LxLy 




XkV-^ 


qi.q2 


LxLy 


VI 


(A4a) 

(A4b) 

(A4c) 


g/[q2 x< qi]z^s/2^(2) 

Then the normalization condition is 

(d'ld') = |d''='‘|2 -q |d'®|2 -q -p |d'(2)|2 = (A5) 

and the system of linear equations can be written as 


/ \ 

Xpe 

X]/^ 

X]/(2) j 


.Ik =0, (A6) 

where the effective Hamiltonian matrix is given by 


H. 


4x4 

eff 


2A„ Fk Vjf 0 \ 

Vk En + E'^ 0 Vk \ 

Ff 0 E^+E'^ Vk ]■ ^ 

0 Fk Vk 2E'^ ) 
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The upper left matrix element of {Vl — ^ coincides 

with Eq. (20), which proves their equivalence. 


For the coupling to T point phonons, we can look for 
the eigenstate in the form 


1^) = q, cr) + ^ q, cr) + 

q,(7 q,cr 

^ 'y ^ ^qi ,cri ,q2 ,CT2 Iftl 5 5 ^^2) 5 

(A8) 


qi,CTl,q2,o-2 


where 'I'qi,cri,q 2 ,o '2 = ^q 2 .o- 2 ,qi,o-i■ The basis states are 
defined as 


h) = /Y ci,p,KC-n,p,K\0), 


|e,q,a) = 


p 

I 2ttI% y. 


|/i,q,cr) = 






^ ^n',p-q,„/2,iC^-",P+9x/2,if^q,(T|0)j 

|qi,cri,q2,cr2) = 


2ttI 


LxLy 


M- ^ ^ ^WQy^B 


(A9a) 


(A9b) 


(A9c) 


^ ^n',p-Q,„/2,if^-’T’'.P+Qx/2,K^qi,cri^q2.o-2l®)- 

(A9d) 


Again, in the last line Q = qi + q2, and we focuse on 
the sector where the initial electron-hole pair \eh) is in 
the K valley. Note that |qi, cti, q 2 , ( 72 ) and |q 2 , (72, qi, (7i) 
represent the same state, so the wave function in the two- 
phonon sector is assumed to be symmetric, g-, .q 2 .o -2 = 
and the factor 1 /2 is necessary in Eq. (A 81 . 




q2,'5’2,qi,'7i 5 


The Schrodinger equation results in the 

following system of linear equations for the coefficients in 


Eq. (A 8 I: 




Xrv'^ 

= -Y^-n',-n{(l) 


q,^T 


YJL'i<DK,<r, (AlOa) 


q,fT 


q,cr 


LxLy 


Xtv^ 


{n-En-En'-ujT)^X,= 


q',<T' 


I-'x^y 


{n-Er,- - ccr)4'E = 


Xrv^ 

= - 


q',<T' 


(fi _ - 2a;r)«' 


qi,^7’i,q2,o’2 




(AlOb) 


(AlOc) 


(AlOd) 


This system differs from Eqs. ( |A3a )-( A4c) by the pres¬ 
ence of two extra terms in the last equation, due to the 
symmetry of the wave function in the two-phonon sec¬ 
tor. These terms correspond to the exchange of phonons 
between the electron and the hole, and because of them 
the whole system becomes non-separable. 


Appendix B: Some properties of Eq. (32) 


In the dimensionless variables z = {e — ojk)/Vk, g = 
VkGu, S = {En—UJK)/VK, Eq. (321 acquires the following 
form: 


-=z-6-^ [ re 
g ^g J zjg - 2 ^”e-«/n! 

0 

The function 


XF n('a) — 


^”e ^/n! 


u - 2 ^"e-?/n! -h i 0 + 




(Bl) 


(B2) 


is real unless 0 < m < = 2n"’e“"'/n!. It behaves as 

l/u at |m| — >■ 00 . At M — >• —O''" it diverges as — ^+(|m|), the 
larger of the solutions of the equation |u| = 2 ^”e“^/n!. 
At w —>■ it has a square-root singularity, provided 

that n > 0. For n = 0,it can be evaluated explicitly: 
Eo(u) = ( 1 / 2 ) ln[it /(u - 2 )]. 
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Treating u = z/g as a new unknown variable, we 
rewrite Eq. (B1) as 


V2 


Fn(u) z(z -5)-' 


(B3) 


so that at any fixed n, it contains only one parameter 


w = z{z—6). Numerical solution of Eq. (B3) is facilitated 
by noting the following properties. 


• At large positive w, correspon ding to large \z\ or 
large negative detuning 6, Eq. (B3) has a real solu¬ 
tion u « z{z—S), corresponding to the unperturbed 
Green’s function g{z) = l/{z — 5). Another real so¬ 
lution is at M « 


• As w is decreased below some critical value , 
the two real solutions merge, and two complex so¬ 
lutions appear [the sign of Imit should be chosen 
to give the correct analytical properties of g{z)]. 

• For w large and negative (i. e., moderate \z\ and 
large positive 6), there is also a real solution u « 
z{z — S) < 0, as well as another one close to u = 0. 


The two solutions merge when w = It is 

easy to see that for w = 0 real solution exist, so 

> 0 . 

• In the region w™” < z{z — 6) < the solution 

is complex, which means that there are two discon¬ 
nected intervals of z where the spectral density is 
non-zero. 

Specihcally, for n = 2 we have w™®" = 0.0296..., 
w^ax — 1.8265..., the corresponding values of u are 
u = —0.070... and u = 0.901.... For n = 0, ic™™ and 
w^‘>^x solutions of the equation ln2u> = 2w — 5, 

= 0.00339 ..., wS®'®’®’ = 3.4684 ...; the correspond¬ 
ing values of u are u = 4wo/(2r(;o — 1) which gives 
—0.0137... and 2.337.... The numerical smallness of 
yjmin gjygg j.jgg sharp but finite features in the spectral 
density. For example, the shape of the spike in Fig. 
is given by 


TT ^z-g{z) 


oc 




Wn 


-z6 


(B4) 


The typical width of the spike is ^ ry™®"/5, which ap¬ 
pears very narrow on the scale of the figure. 
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